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Abstract. In the present paper the macroscopic limits of the kinetic model for inter- 
acting entities (individuals, organisms, cells) proposed in [3] are studied. The kinetic 
model is one-dimensional and entities are characterized by their position and orienta- 
tion (+/-) with swarming interaction controlled by the sensitivity parameter 7. The 
macroscopic limits of the model are considered for solutions close either to the dif- 
fusive (isotropic) or to the aligned (swarming) equilibrium states for various 7. In 
the former case the classical linear diffusion equation results whereas in the latter a 
traveling wave solution does both in the zeroth ('Euler') and first ('Navier-Stokes') 
order of approximation. 



1 Introduction 

The paper [5] focused on biological alignment (swarming) characterized as co- 
herent motion of groups of entities (individuals, organisms or cells) into the same 
direction. We understand alignment as an adaptation of orientation to that of 
neighborhood. It results from the interplay of the behaviors of entities by local 
interactions. Swarms, also called herds, flocks, schools, clusters depending on 
whether they refer, respectively, to insects, mammals, birds, fish, bacteria and 
cells are often observed in nature. The typical examples are herds of sheep, 
flocks of birds or schools of fish (see References in [3] and [11] [191 HSl [30]). In 
some of these phenomena alignment is a typical behavior. 
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The phenomena of ahgnment (swarming), viewed as a cooperative phe- 
nomenon arising from the interaction of a number of entities, can be analyzed by 
mathematical models. The entities are characterized by their position, orienta- 
tion (velocity) and migrate in space and time. Various microscopic, mesoscopic 
and macroscopic models describing swarming have been proposed. The dynam- 
ics of alignment in the space homogeneous case was considered e.g. in [T^[3T]. 
Microscopic models, for example the cellular automaton approach, [11 [11], and 
a simulation model by Ben-Jacob et al. [S| allows to distinguish individual 
entities. By means of a macroscopic model based e.g. on hydrodynamics (see 
|5S]), it is possible to study the dynamics of densities of entities. The review 
of various individual-based models based on ODEs is presented in [5]. In a 
number of papers the authors have studied stochastic models for swarm ag- 
gregation, where the individuals, in addition to the classical Brownian random 
dispersal, are subject to long range attraction and short range repulsion. For 
instance, see [S] and references therein, under suitable laws of large numbers it 
was shown that, for a large number of individuals, the evolution of the empirical 
distribution of the population can be approximated using a suitable limit non- 
linear degenerate and nonlocal parabolic equation. In [M] a class of stochastic 
individual-based models, written in terms of coupled velocity jump processes, 
was presented and analyzed. Moreover, in the limit of large populations, a sys- 
tem of two kinetic equations with nonlocal and nonlinear right hand side was 
derived and analyzed. 

Mesoscopic models have been considered e.g. in [18l[3l[24l|9l[T]; see references 
therein. Such models are usually of Boltzmann type, i.e., they are related to a 
statistical description of one test entity. 

The relations between the microscopic, mesoscopic and macroscopic mod- 
els is one of the most challenging mathematical problem. Its importance is 
particularly visible in justification of the macroscopic models by a well-defined 
behavior of the entities of the system - see e.g. p[7ll8lfT0 l [T4 l [T5 l fT6 l [20 l [22] 
and references therein. 

In the present paper we study the macroscopic limits for the discrete ve- 
locity Boltzmann-type (mesoscopic) model based on a well-defined microscopic 
'majority-choice' interaction introduced in [5]. We are interested in the forma- 
tion of alignment of the population of entities that do not possess a leader but 
polarization arises as a result of local alignment interactions. 

The model, proposed in [3], is one-dimensional and the entities are charac- 
terized by their position and orientation (-1-/-) while the alignment (swarming) 
interaction is characterized by a sensitivity parameter 7. In the present paper 
the macroscopic limits of the model are considered for solutions either close to 
the diffusive (isotropic), or to the aligned (swarming), equilibrium states for 
various values of 7. In the former case, the classical linear diffusion equation 
results whereas in the latter we obtain a traveling wave solution both in the 
zeroth ('Euler') and the first ('Navier-Stokes') order of approximation. An in- 
teresting feature of the considered model is that the nonlinear kinetic model in 
the 'hydrodynamic' limits gives, in both cases, linear macroscopic equations. 
Moreover, it follows, at least formally, that all terms of the asymptotic expan- 
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sion of the bulk part of the approximation, which are complementary to the 
traveling wave solution (its' 'kinetic^ part), vanish and thus the bulk part of the 
approximation only consists of the traveling wave, whereas the complementary 
part is made up only of the initial layer terms. 

The plan of the paper is as follows. Section [5] presents the model that was 
introduced in [3] . In Section [3] wc prove a unified existence and uniqueness 
result valid both in weak and strong interaction cases. Section |4] deals with 
formal macroscopic limits of the model. This section also indicates some inter- 
esting open problems: the hydrodynamic limit in the 'diffusive' case. Finally 
we propose an asymptotic result in the 'aligned' case. The result is proved in 
the spirit of the Tikhonov-Vasileva theory (cf. [57]) that was formulated only 
in the finite-dimensional dynamical system case. Therefore it is an example of 
a generalization of the theory available for ODEs into some infinite dimensional 
case. 



2 The model 

Let / — f(t,j,x) be the probability that an individual is at time t at point x 
with orientation j, where i > 0, a; G R and j E { — 1, !}• Here the orientation 
plays the role of discrete velocity in the sense that it give the direction of motion. 

We describe the migration together with the changes of the orientation vari- 
able where the latter is modeled by an interaction operator. We have the fol- 
lowing general expression (valid for every migration- interaction case): 

f{t + At, J, X + jAx) - f{t, -J, x) ^{t, -J, x) + fit, j, x) q3'(t, j, x) , (1) 

where 

*P(f, j, x) = Probability ^a change of orientation in At (t, j, x)^ , 

is the probability of a change of orientation in the interval of time [t, t 4- At] of 
an individual that is at t at position x with orientation j and 



'9'{t,j,x) = 

Probability ^no changes of orientation in At 
l-^{t,],x). 



Different choices of the probabilities *p give rise to different models, see e.g. 
[5] where, in particular, wc proposed 

x( E f{t,k,x + al)>0)( j: f{t,-j,x + al)y 

\ U I Ul / \ / Li / 

W, J, ^) = /' At , (2) 

( E /(i,-J,.T + aO +( E /(i,J,a; + a/) 
^i=±i ' ^i=± ' 
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where a > and 7 > are parameters and ^ /(<, j, x+al) is the neighborhood 

i=±i 

density in direction j, x(true) = 1, x(fa'lse) = 0. The parameter 7 describes 
sensitivity of the interaction. If 7 is smah (close to 0), then the probabihty 
of a change of orientation only weakly depends on the actual orientation. On 
the other hand, for large 7 the probability of a change of orientation strongly 
depends on the actual orientation. 

Assuming that a = Ax and Ax = j At, in the limit At — >■ 0, we obtain the 
following system of two equations 

x{ E fit,k,x)>o) 

dtf{t, j, x) + jd^f(t, j, x) '"'^'^ 



{.f{t,-3,x)y + {f[t,3,x)y 

X (^f{t,-],x)i^f{t,j,x)y ~ f{t,j,x)(^f{t,-j,x)y^ , 3 = -1,1. 

(3) 

For 7=1, Eq. ([3]) decouples into two free-streaming equations. Throughout 
the paper we consider either 7 €]0, 1[ or 7 > 1. 

The non- negative and non-zero equilibrium solutions /, corresponding to the 
space homogeneous version of Eq. ([3]) , can be found from 

/_i(/i)^-/i(/_i)^ = 0. (4) 

The only non-zero solutions / = (/i, /-i) of Q are given by 

/i = /_i>0, (5) 

or 

fj = and f-j > , for some j e {-1, 1}. (6) 

Equation ^ corresponds to equal probabilities of both orientations (a dif- 
fusive picture), whereas Eq. (|6]) is related to an aligned picture. 

In the space-homogeneous case the trajectories of the corresponding ODEs 
are contained in the straight lines defined by 

/-i+/i=c, where c = /_i(0) + /i(0) > , (7) 

and have different types of behavior according to the value of 7. For 7 > 1 

lim /- 1 = , lim A = /i (0) + /_ 1 (0) , for A (0) > /_ 1 (0) , (8) 

lim /_i =/i(0) + /-i(0), lim/i=0, for/i(0)</-i(0), (9) 

t—>oc t^oo 

whereas for < 7 < 1 we have 

, , /i(0) + /-i(0) 
hm /i = hm /_i = . (10) 

t—>oo t^OO Z 
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The important problem whether the similar behavior can be observed in the 
spatially nonhomogeneous case remains open. 

We note that the spatially nonhomogeneous Eq. ^ admits simple solutions 

f{t,j,x) = f{t,-j,x) = const. > 0, (11) 

and 

fit,j,x) = 0, fit,~j,x)^^{x+jt), J = 1,-1, (12) 

where (j) is a given non-negative function, obtained by equating both sides of 
([3]) to zero. 

Thus we observe that the possible asymptotic behavior for Eq. (jSj may differ 
from that of the Carleman-type equations studied in [531 1131 IHl US] • 



3 Existence 

The (global) existence and uniqueness results for 7 > 1 as well as the existence 
of an entropy functional were established in [3] taking advantage of the Lipschitz 
continuity and the conservativeness of the system. 

Here we propose the global existence and uniqueness results in an appropri- 
ate Banach space for both 7 > 1 and 7 < 1. We will consider one dimensional 
periodic boundary conditions, that is, we assume x G T, but similar results are 
possible for x g T'', where d> 1 and T"* is a o?-dimensional torus. 

Let Xi be the space of real-valued functions equipped with the norm 



/lli= E / 

^ Li 



The cone of non-negative functions in Xi is denoted by X^ . We define the 
following operators 

Q[f] = Q^[f]-Q-[f], 

where 

A:=±l 

and 



fc=±i 



Note that, for / > 0, 
where 



P{j,x)+p{~j,x) ' 

/(j,x)r(-j- x) 

Pij,x) + Pi~j,x) • 

Q[f]=R[.f]-P (13) 



i?[/](.-,x)=r(,-.)-^(^'^)+^(-^'^-) 



Pil,x) + pi^l,x) 
We denote 

h'^ij,x)h{-j,x) 



,5)(j>') 



h^ij,x)+g^{~j,x) ' 
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and 



Let 

fHt,j,x) = fit,j,x + jt). (14) 
We define the following sequences {gn}neN and {/injneN: 

50 = 0, 

while ho is the solution of the problem 

hl{t,j,x) = hl{t,-j,x), j = ho\t=o=fo, 

where /o € is the initial datum. 
Moreover, 

gl+gl = 0, gi|t=o = /o, 
+P«(/i„_i,5„_i)5«, = Q«(5«-i,/in-i) , g„|t=o = /o , 
for n > 2, and 

/il + P«(5n-1, = Q«(/l„_l,5„-l) , /i„|t=o = /o , 

for n > 1. 

Therefore we have 

gi{t,j,x) = fo{x- jt)e~* 
\\ho{t)\\i < ll/ollie*, 

and 

< gi{t,j,x) < fo{x-jt) < hi{t,j,x) < ho{t,j,x) . 
Assuming that 

fo{x) > M a. a. X G T (15) 

where /i is a positive (> 0) number, we obtain that on any compact interval 
[0,T] 

< /ie ^ < gi < g2 < 9z < ■ ■ ■ < h-i < hi < hf) . 
Therefore, by the theorem of Bcppo-Levy both sequences converge in Xi 

9nt g, hnlh, 

in Xi and 

9<h. 

On any fixed time interval [0,r] wc may apply the inequality 

0<h-'-g-' <^ifie-^r-\h~g) 

to show that g = h. 



Thus f -.^ g = h is a mild solution to the Eq. © on [0,T], T > 0, in the 
following sense 

/eL-(0,T;Xi), 
t t 
f\t) + / /) fl = /o + / Q«(/, /) ■ 



for a. a. t G [OiT]- The similar argument as the one used for g ^ h shows that 
the solution is unique. Moreover, the solution is in . 
Therefore we have 

Theorem 3.1. Let fo G be such that 1115]) is satisfied for some positive fi. 
Then, for any given T > there exists a unique mild solution f ~ f{t) (z Xi of 
Eq. (0) with the initial datum /o on [0,r]. Moreover, 

• f{t) e X+ for a.a. te [0,T], 

• t^e-^ < f{t) for a.a. t G [0,r], 

• ||/(t)i|i<||/o||ie^/ora.a. t e [0,T]. 

Remark 3.2. Note that if ([T5|) is satisfied, then / is a mild solution of Eq. ^ 
in which each term is on the right hand side is separated from zero and thus 
Frechet difFerentiable. Hence, the solution becomes a classical solution provided 
the initial datum is differentiable. 



4 Formal macroscopic analysis 

In this section we follow the strategy developed in the case of the so-called 
hydrodynamic limits for kinetic equations in [T5| . 

It is well known that the macroscopic limit is defined by e — >■ 0, where e > 
is a suitable small parameter (the Knudsen number). Classically, see e.g. [T5] . 
Equation ([3]) can be considered in the following two dimensionless forms 

[dt+jd,^fit,j,x) = ^Q[f]{t,j,x), (16) 

and 

(edt+jd,)f{t,j,x) = ^Q[f]{t,j,x) , (17) 

where the nonlinear operator Q is given by the RHS of Eq. ((S]) . Equations (|16p 
and ([T7| with the initial data 

/ (18) 
lead to singularly perturbed problems in the limit e — > 0. 
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We introduce some basic notation. The 0-th order term in the Hilbert 
expansion is the function M (a Maxwelhan) which satisfies 

Q{M) = 0. (19) 

Therefore 

Mit, j, x) = T], g{t, x), ryi + r/_i = 1 , (20) 

where either 

Vi^V-i = l, (21) 

or 

77fc = 1 , ry-fc - , (22) 

for some k g {—1,1}, and g is the macroscopic density. Therefore we conclude 
that two macroscopic hmits can be considered for initial data close to (PO)) either 
with ((2T|) or (|22|) . This is, however, strictly related to stability properties of the 
solutions. 

Consider the first Frechet derivative (linearization) L of the operator Q at 
M. It has the form 

1-7 



if (EI]) holds, and 



LfU) = ^{.f{-j)-f{j)) (23) 



\-.fi-k) forj = -fc ^ ^ 

if p2|) is satisfied. For simplicity the dependence of t and x is not explicitly 
indicated. 

The null-space of L, called the hydrodynamic space, is spanned by the Max- 
welhan M. Since the aim of asymptotic analysis is to isolate the slow dynamics 
in the hydrodynamic space from fast, transient, behavior in the complementary 
space, the best way is split the phase space into subspaces reducing L. Hence, 
we call the kinetic space the space spanned by all remaining eigenvectors of L 
(corresponding to eigenvalues with negative real parts). The kinetic space will be 
denoted by W. This can be done by spectral projections. For a general matrix A 
having as the dominant simple eigenvalue with the corresponding eigenvector 
with e and e* being the eigenvector of the transpose A* corresponding to and 
normalized so as e • e* = 1, the spectral projection V onto the null-space of A 
is given by 

Vy= (e* ■y)e. 

Then clearly the complementary projection is defined as = I — V , where 
/ is the identity matrix. Specifying this theory to the case in hand, we see 
that in the diffusive case ([23]) the null spaces of L and L* = L are spanned 
by, respectively, (1/2,1/2) and (1,1) (where we have taken into account the 
normalization). Thus, in this case we have 

P(/a),/,-i)).fM±im,M±Zhl)) (.5) 
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and 

(/(I). /(-!)) ^(M^.itM). M 



The situation is, however, different in both cases described by ([22|) . For fc = 1 
the eigenvector of L corresponding to is given by (1,0), whereas (0, 1) is the 
respective eigenvector for fc = — 1. For the transpose matrix the eigenvector is 
given by (1, 1) in cither case. Thus, we obtain 



and 

if A; = 1, and 
with 



7'(/(l),/(-l)) = (/(l) + /(-l),0) (27) 
P^(/(l),/(-l)) = (-/(-l),/(-l)), (28) 

= (0, /(l) + /(-l)) (29) 

7'^(/(l),/(-l)) = (/(l),-/(l)) (30) 

for k ~ —1. We observe that the above expression give a prehminary justification 
for the terms diffusive and ahgned picture. Indeed, in the first case the total 
density g = /(l)+/(— 1) is uniformly spread between the two available directions 
of motion whereas in the other two the total density is concentrated along one 
of the directions of motion. Let us introduce an operator D by the formula 

D{f{l, x), /(-I, x)) = (5t/(l, x) + d^fil, x),dtf{-l, x) - dj{-l, x)) . 

Referring to the problem we note (see [TS]) that the abstract zcroth order 
('Euler') macroscopic approximation is given by 

VDM^O, (31) 

whereas the first ('Navier-Stokes') order macroscopic approximation is given by 

VDM = -eVDV^L-^V^DM (32) 

in terms of the macroscopic density g. 

Specifying for our three cases, we obtain the following formulae 

{^{dtg + dxg,dtg — dxg) in the diffusive case 
{dtg + dxg, 0) in the aligned case, k = 1 

(0, dtg — dxg) in the aligned case, k = — 1, 

which, upon application of an appropriate V , yield the following Euler approx- 
imating equations for the diffusive and aligned regimes with k ~ 1,-1, respec- 
tively, 

dtg = , 

dtg + dxg -0, (33) 
dtg- dxg =0. 
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To derive the Navier-Stokes picture, first we find 

{{dxQT—dxO) in the diffusive case, 
(0, 0) in the aligned case, k = 1, 

(0,0) in the aligned case, k = — 1, 

Clearly, the operator L is not invertible but becomes such if reduced to the 
kinetic space W which, in the diffusive case, is spanned by (1, —1). We see that 
V^DM e W and the solution to L{y,z) = {dxQ, —dxg) subject to y = —z is 
{-dxg,dxg)/{l-'-f)- Next, 

DV^L-'V^DM = (1 - j)-\-dtdxg - dig, dtdxg - dig) 

and finally 

VDV^L'^V^DM = -{l--/)-\dlg,dlg) 

which, combined with VDM = {dtg,dtg), yields the first order macroscopic 
approximation (|32p in the form 

dtg==-^dlg. (34) 

1-7 

Clearly, the Navier-Stokes approximation in either aligned case is again the 
respective Eulcr approximation. 

Next we consider the diffusive picture regime for problem the parabolic scal- 
ing (jl7|) . Performing similar considerations we see that the zeroth order macro- 
scopic approximation is given by 

dtg^j^dlg. (35) 

1-7 

Equations ([M)) and ([55)1 are linear diffusion equations for < 7 < 1, with the 
diffusion coefficients and , respectively. Linear dependence on e of the 
diffusion coefficient in the limit equation ([M)) in the case of scaling (ITS]) is typical 
(cf. [Ullin]). In the case of 7 > 1 Eqs. and are the (linear) backward 
diffusion equations. The backward diffusion equation is less mathematically 
tractable than the diffusion equation — see however [T71 Example 8.6]. 

Thus, we cannot expect any stability properties in the diffusive regime if 
7 > 1. On the other hand, one may expect that for < 7 < 1 the asymptotic 
relationships between the solutions of ()16p or (fT7|) and (p4| or (|35|) , respectively, 
hold. 

Next, let us consider the aligned regime for Eq. with initial data ([T5)) . 
Let k € {—1,1} be fixed and, to avoid technical difficulties we consider only 
7 > 1. We have seen that the macroscopic limit at both Euler and Navier- 
Stokes levels does not depend on the choice of 7 and is given by 

dtg + kdxg = 0. (36) 

Equation p6p is a linear first order equation and for any smooth initial datum 
F its unique solution is given by 

g{t,x) ^ F{k,x ~ kt) , t>0, x€T. (37) 
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Clearly, the solution inherits all properties of F such as boundedness, differen- 
tiability, etc. This solution is a traveling wave (in the direction k). 

Next we investigate in more detail to aligned picture and, in particular, we 
provide a more refined expansion formulae. To simplify notation, let us fix 
fc = 1; the calculation for fc = — 1 being symmetric. 

Then, following ([?f)) and ([25)) . we have the hydrodynamic part of the solution 
given by p = /i + /_i and the kinetic part defined as w = /_i so that 

(/i,/_i) = (l,0)e+(-l,lK (38) 

Using fi = Q — w and /_i = w, we see that the system pB]) can be written in 
hydrodynamic and kinetic variables as 

dtQ + dxQ - 2dxW 
dtw — dxW 

First we look for the bulk approximation. 

We perform the asymptotic expansion in the spirit of the Chapman-Enskog, 
that is, we do not expand the hydrodynamic part and seek approximation 
[g, w) — {g, wq + ewi), see e.g. [20]. This produces 

dtg + d^g - 2dx{wo + ewi) = , 

edt{wo + ewi) - edx{wo + ewi) (40) 

_ (g - [wq + £Wi))(wo + £Wiy - [wq + £Wi)(g - (wq + £Wi))'^ 

(g - {wo + £Wi)y + (wo + swiY 

Taking 0(1) terms on both sides, we obtain 

Q _ (g - wq)w^ - wo{g - wq)'^ 
(g - wo)'y + 

which yields wq = or wq = g. If we recognize that the kinetic part for k ~ 1 
is the hydrodynamic part for k = —I, then the latter option simply give the 
fc = — 1 case and we are left with wq — 0. To get wi we write 

(g-eWl)(£U;i)T -£Wl(g-£Wl)T _ -£g''Wi + 0(£"""^2,7}) 

(g-ewi)^ + (ewi)^ ^ gT(l - e-fwig-^ + 0(£™i2.7})) 

= -eWi + 0(£"""{2,7}) 

and, since the remaining LHS terms in ([40| are 0(£^), it follows that wi — 0. 
Hence, we recovered the relevant, second, equation from p3[) . Wc note that in 
this model and the presented level of expansion Chapman-Enskog and Hilbert 
methods yield the same approximation. 

Remark 4.1. By induction wc can prove that all terms of the expansion w = 
Wo + £Wi + . . . + £""Wn + ■ ■ . are zero. Indeed, if Wk = for fc = 0, 1, . . . , 7i — 1, 



0, 

1 (g — w)'w'^ — w{g — w)'^ 
e {g-w)-i+w^ 



(39) 
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the calculation as above gives 

(g - £"W„)(£"W„)T - e"-Wn{Q - g"Wn)^ 
(£--£"W„)T + (e»W„)7 



F^o'ly, _|_ nCpmin{2n,7n}\ 



£>T(1 - e"7W;„£»-l + C)(emin{2n,7n})) 

and since the LHS of PU)) now is 0(£"+^), we find !/;„ = 0. 

Next we incorporate the initial layer correction and look at the approxima- 
tion {g{t),'w{t)) [g[t) +r(-), h{^)), where g is the solution to the initial value 
problem 

dtg + 0^9 = 0, giO,x) = go{x)^F{l,x). (41) 

Then, denoting r = -, the initial layer corrections r, h are obtained by an 
asymptotic expansion of 

dtg + e-^drr + d^{g + r)^2d^h = 0, 

{g + r - h)-* + hf ^ ' 

Since we expect r to decay exponentially to as r — > cx), from the first equation 
we obtain r = 0. Then, setting £ = in the second equation and taking into 
account that g{t^x)\e=o = g(eT)|e=o ~ Qo we obtain 

[go - hp + 

h{0,x) = F{-l,x). (44) 



5 Aligned picture 

We consider the 'aligned picture^ for Eq. in the case 7 > 1. Wc assume 

that k G {—1, 1} is fixed. 

Let X G T, where T is the 1-dimcnsional torus. All functions are then 
interpreted as periodic functions on the unit interval [0, 1]. 

Let Xoo be the Banach space of continuous functions defined on T with the 
norm 

ll/lloo = • 

je{-i,i} 

We consider Eq. i.e. 

s(dtf,+jd.,f,) = "^^^'^[-'^ °^ (f-jf] ~ f,fl,) , J^k,-k, (45) 

j ■'-3 

where fj{t,x) = f{t,j,x). 
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Following the formal results derived in the previous section, we are looking 
for the solution to Eq. p5)) in the following form 



fk(t,x) = g{t,x)-h[-,x] +eu{t,x) , f-k{t,x) = h ( - ,x ] + ev{t,x) , (46) 



where (ew, ev) is the error of the approximation. We recall that g = g{t, x) is 
the 'bulk solution' , h = h{T,x) is the initial layer solution' term, r = | is the 
stretched time variable. Then the error satisfies 

e^iptu + kdxu) — drh — ekdxh 

x{g + eiu + v)>0) (^{h + ev) {g - h + euy - [g - h + eu) [h + evy'j 
{g - h + euy + {h + evy 
[dtv — k dxv) + drh — ekdxh 

x{g + e{u + v)>0)(^{g- h + eu)(h + evy - {h + ev) {g - h + euj'^'j 
[g — h + eu)^ + (h + ev)'^ 

(47) 

where g satisfies ([37)) . h — h{T,x) satisfies pS]). p4| and go{x) = g{0,x) = 
F{k,x). 

First we consider the initial layer problem for h. We have 
Lemma 5.1. Let k G {—1,1} and the initial data F £ X be such that 

F{-k,x)>fi, VxeT, (48) 
for some /i > 0. There exists c-y > 1 such that, if 

F{k,x) > c^F{-k,x) Va-eT, (49) 



then there exists a unique solution h — h{T,x) to the problem ^43y, !i44\ ) for 
any t > and 

<h{T,x) < F{-k,x)e'^^ Vr>0 VxeT, (50) 

for some 6 > 0. Moreover if, additionally, F is continuously differentiable and 
Cy is sufficiently large, then the solution h = h{T,x) to \4S^ , continuously 
differentiable and 

\dxh{T,x)\<const.(\\F\\^ + \\dxF\\^y~^^^ Vr>0 V;reT, (51) 
for some 5i > Q. 
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Proof. Problem (|33]), (pi)) is an ODE problem, where the variable x is a param- 
eter. For 7 > 1 the operator defined by the RHS of (|43)) is Lipschitz continuous 
on C(T). Under the condition dH]) for a.a. a; e T and r = the RHS of Eq. 

is negative. Therefore h strictly decreases with respect to the variable r 
for all X G T. Now, 



■(^{go-h)h^ -h{go-hy') 



(go - hp + /iT 

[go - h)-' + 

Since, by 

h F(-k) 1 
< i < , 

ga-h - F{k)- F{-k) - Cy-V 

wc can write 

(go - M((go - hp-' - h''-') ^ 1 - ^^-^ 
{go ~h)i + h'r 1 + t?T 



d>0 (52) 



uniformly in x € [0,1] and r > 0. We note that 6 can be made as close to 1 
as we wish by making sufficiently large but since close to the equilibrium 
h = the left hand side of ([5^ is close to 1, (5 could not exceed 1. This leads 
to the conclusion that h decays exponentially to for all a; G T and the rate 
of convergence can be controlled by e~^'^ uniformly in x G T. Moreover it is 
strictly positive. Thus the unique solution of Eq. p3)) satisfies ((50| . 

It is evident that if F is continuously diffcrentiable with respect to x, then 
also h is. Wc denote h' — d^h and £Iq = dxgvi- The function /i'(t, x) satisfies 

p.,, l{Qo-h)'^-^{g'o-h') + -ih'i-^h' / ^ ^ 

drh^ ((g, _ + fe7)2 {9.-h)h[[go-hr -h? 

^ -//(go - fe)((go - - /l^-l) - fe(g[, - feO((go - h)-^-^ - 

(go - h)i + /iT 
^ ^ fe(go - hr-\do ~ h') - h-'-^go - h)h' 

(go - hp + h-y 

The first term on the RHS of ([55)1 docs not cause any difficulty: the condition 
(|49|) is sufficient to ensure the negativity of the term multiplying h' . Here we 
are using ([SO]) . Considering the second and the third terms on the RHS we see 
that if 

-(go - 2h) ((go - hp-' - /i^-i) + (7 - l)higo - hp-' 

+ {^-l)h^-\go-h) <0 (54) 



for all X G T, then ([5T|) is satisfied. It is easy to see that if c-y is sufHciently large 
(for example C2 > 6/(3 — -y/S)), then ([54|) follows. This ends the proof. □ 
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We may state the following result 



Theorem 5.2. Let the initial data F G X^o be nonnegative functions with 
continuous second derivatives and such that |^<$[ ), for some fi > 0, is satisfied. 
For any T > there exists Eq > and > such that if 



min F(k,x) > c^ niaxi<"(— fc, a;) 



(55) 



is satisfied for given k g {—1,1}, then for e € [0,eo[ the Cauchy Problem 
for Eq. ( |-^5[ j, with the initial datum F, has a mild solution f = f{t) in Xoo on 
[0,T]. Moreover, 



sup 
te[o,T] 



fit) - Qit) 



< ctE . 



(56) 



where 



h = h 



Q{t,k,x) 
Q{t,-k,x) 

h(|,-fc,.T 

j = k,—k, 

h{0, x) = F{^k,x) VxeT. 



F{k,x~ kt) , 
0, 

-h[lx), 



Proof. Let T > be fixed. The proof follows by the analysis of the system (|47p 
that, by P5)) . is equivalent to Eq. The methods of Theorem 13.11 may be 

used to show the existence and uniqueness of the mild solution / = f{t) in X^o 
for any fixed e > 0. The solution satisfies 

• f{t) e Xryo and f{t) > for a.a. t e [0,T], 

• A^e"? < fit) for a.a. t e [0,r], 

• ll./(t)lloo < !l./o!|ooe? for a.a. t e [0,T]. 

Wc may rewrite Eq. (j45p in the following form 

( ^^^+^;^^^^-^^ (i7'^-/27^)y, .i-k,-k, (57) 

where " jj" is given by ((Ti| . Under the assumption the functions fit,k, x+kt) 
and /(t, —k,x — kt) are increasing and decreasing functions of t (for any fixed 
x), respectively, and we have 

fit,k,x) > F{k,x-kt) > minF{k,y) > 

yeT (58) 
> Cj max F{—k,y) > c-yF{—k,x + kt) > c-yfit,—k,x), 
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for any t > and any x € T. 

Now we assume (|46|) and consider Eq. (|47)) . taking in consideration that 
^ + eu + ef > 0, we obtain 

dfU + k dxU = drh H — kdxh 

(h + ev) {g — h + eu^ ~ (^g — h + eu) (h + evY 
e"^ [{g - h + euf + {h + evY^ 

dtv~kdxV = — t: drh -\ — kdxh 

e 

^ [g — h + euji^h + evY — [h + ev) {^g — h + euY ^^^^ 



(^{g-h + euy + [h + evy'j 



Therefore, using (03]) and P^, 

dfU + k dxU = 



1,„ , 1 / (fi - + £«)'''(£» + eu + eu) ig-hYg 
-kdxh + ^ ' 



e^ y /i + eu)^ + (/i + eu) [g - hy + h^ ) 

1 1 ;i(go-/i)''-((?o-/i)/i''^ 



+ {ga-hy + h-y 



dfV — k dxV 



1 if {h + evy {g + eu + ev) K^g 

—kOxh H ^7 



V 



{g~ h + euy + {h + ev) [g - h)^ + h^y J 
1 1 ( {g- h)hy - h{g- h)^ {go - h)h'' - h{go - h)' 



'V {g-h)^ + hf {go~h)^ + hf J 

(60) 

Assume now that g, h, u, v are given and consider the foUowing functions 

„ {g - h + eu)'^ {g + eu + ev) 

^'^'^ ^ {g-h + euy + {h + evy ' 

and 



S2(e) 



{h + ev)^{g + eu + ev ) 



[g — h + eu) + (ft, + ev) 
We have 

1 

S,(e)-S,(0) = e J E'^{9e)de, i = l,2. 



The derivatives S^(e), i = 1,2, may be written as the sum 
s;.(e) =r,4(e)u + r,,2(e)u, 
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where, denoting Q = g — h + eu and A = h + ev, 



ri,i(£) 
ri,2(e) 
r2,i(e) 
r2,2(e) 



927 _ jeiA''-^{g + eu + ev) + e'>A'> 



2 

07 + At ' 



97 A-/ - jQf-^A^ {g + eu + ev) + A^t 



2 

97 + At ' 



e^AT + 7eTAT-l(^ + £li + £u) + A^T 



2 

97 + At ' 



Let us consider 

1 I A7 e+eit+Ef I A7 

-l-iaiEj - ; ri 

1 + Ai 



, ,A'f-^ g + eu + ev At-i\ / At 

1 + (h + ev)- • ^ ' ' 



97 9 97 ; \ 97 

Now, by dH]), dnSl) and dSS]) we get 

g — h + eu > min{g — h, fk} 

> min{F(fc) - F{~k), F{k)} > F{-k) niin{c^ - 1, c^} > ^i(c-y - 1) , 



and 



(61) 

A h + ev max|F(-fc), /i| 1 , , 

= T < -, ; / < ^ 0, c^^oo. (62) 

e p-h + eu- (c-^-l)F(-fc) - ^ ^ ' 

Performing analogous estimates for other F^j , we see that 

ri,i(e) + r2a(£) = 1, 
ri,2(e) + r2,2(£) = 1, 

|ri,i(£)| < l + ai,i, \y2.i{e)\ < fT2,i, 

|ri,2(£)| < l-CTl,2, |r2,2(e)| < fT2,2, (63) 

where the positive constants cr^j-, i,j = 1, 2, may be as small as we want if is 
sufficiently large and also 

|ri,i(£)| < l + al,{\h\+e\v\), 
|r2,i(£)l < <J2,i{\h\+e\v\), 

|r2.2(e)| < <72A\h\+e\v\), (64) 
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where again the positive constants cri i , 1 1 "'2 2 ^^^^.y be chosen as small as we 
want by taking sufficiently large c^. 

In order to estimate the last term of Eq. (pO]) we assume that r is fixed and 
consider 

^ igjt) - hiT))h^T) - hiT){git) - h{T)r 
^ {g{t)^h{T)r + h-f{T)) 

By (Ha), 

SO that 



{g{t) - /i)T + h-r {go - h)'r + /if ' 
and thus, using differentiability of the initial conditions, ([58| and (|6T|). 



i2 + j)\\F\\l 

//2(c^ - 1) 

Then we have 



< th^\\F'\\^ \J^j'^^ =C^th^. (65) 



dtu + k dxU H — u 

£ 



1 ± 1 

ifc5^/i + J y ri,i(6'e)d6'+ ^ y ri,2(6'£)d6' + ^ j '^'{Ot)dO, 



dtv — kdxV-\ — V 
e 



1 

= —kdxh 

e e 



j r2,i(fe)d0+^ j r2,2(Md^-^ j -^'{0t)A9. (66) 



Integrating along characteristics and using and (|64p . we obtain 
t ^ t 

\\u{t)\\oo < 7/ + y e-t*'||^.(s)|Uds 



t t 

ll^(*)llooll"(*)llaods + — ^^J\\v{s)\\oods, (67) 



t 



t 



' e-'^\\u{s)\\oods + ^le-^\\v{s)\Uds, (68) 
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where C denotes a generic constant independent of e and 5 < min{7(5, 5i\ > 0, 
see Lemma [01 and ((65l) . For convenience of calculations, we fix 5 < 1. 
We can re-write (1551) as 



\Ht)\\^<C + <72.i\\u\\^^t + <y2.2\Ho.,t, i€[0,T], (69) 

where ||u||oo,t = sup ||u(s)||tx)- Taking large enough for 172,2 to be smaller 

se[o,t] 
then 1, we obtain 

IM|oo,t < C(l + c72,i)||w||oo,t, ie[0,T]. (70) 



We may now refine the estimate ((68l) . using (|64)) and (jSOp . with 5 replaced 



t 

C 



ds 







t t 



+ ^je-'^e-i'\\u{s)\Uds+^Je~'^e~i'\\v{s)\Uds 


t t 

+4i/e"'^lk(s)lloo||«(s)|Uds + ai.Je-^|l«(s)||^ds. (71) 





Therefore we find 

t 



- I \\v{s)\\cx.ds < C + 0-2,1 ll^lloo,* + 0-2 2lkl|oo,t 



+ <i(r + e)\\u\Ut\\v\Ut + 42(r + e)||«||L,f (72) 
By dlTD, (Uni) and dZH) we obtain 

llu(i)lloo < C+^\\u\\^,t + ^[iTC\\u\\oo,ti^ + a2.i\\u\\o,.t) 
d 

+ (1 - 0-1,2) (C + 0-2 i||u||oo,t + C'o-j 2(1 + a'2,l||u||oo,t) 

+0^2 i(r + e)C||u||oo,t(l + 

+ 42(r + e)C2(l + a2,i||u||oo,t)') 

so that 

\\u\Ut<C + a^''>\\u\Ut + a^'^\\u\\l^f (73) 
for any t e [0,T]. where o-'^-* and o-*^^' are given by 



cr^*' = + o-^ ;^rC+ (1 - o-i,2)(a-2 1 -f Co-2,10'2.2) 



(1) 

+ C(r + £o)«l+2C42'^2,l), 
= (Ti^ifT2,irC+ ((7^,10-2,1 +Co-^^24l)C(T' + '^0). 
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and tend to with (Jij- — > and cr^ ^ — )■ (that is, for c-y — > oo). Thanks to this 
wc may choose large enough for cr'^^ < 1 to obtain 

\\u\Ut<Cil + a^'^\\u\\l^,), (74) 

for any t E [0, T], where C is a constant. 

Let y := ||u||oo,t, a := C and /? := Ca'-^'K Then, by ([711), we have 

0<-y + a + l3y^. (75) 

Considering the function h{y) = a — y + Py'^, we may choose a > and 
/3 > (a/3 < i) such that 

h{y min) < , for y,nin = ^ ■ (76) 

Moreover = and y cannot experience jumps. Therefore, in particular, 
we have 

y<^, (77) 
and ||w||oo,T is bounded. Then ([7n)) ends the proof. □ 
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